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SECTION  1 
INTRODUCTION 

By  the  very  nature  of  their  reason  for  existence, 
military  aircraft  are  susceptible  to  damage.  The  performance 
of  a  damaged  aircraft  can  be  affected  in  a  variety  of  ways 
depending  on  which  of  the  aircraft  systems  are  affected  and  to 
what  degree  the  damage  is  inflicted.  Structural  damage,  in 
particular,  has  been  responsible  for  the  degradation  of  the 
performance  of  many  aircraft,  and  in  a  significant  number  of 
cases,  for  the  loss  of  aircraft. 

It  is  apparent  that  structural  vulnerability  is  a  highly 
important  consideration  in  assessing  the  overall  vulnerability 
of  aircraft.  One  way  of  assessing  the  structural  survivability/ 
vulnerability  characteristics  of  aircraft  is  to  perform 
sophisticated  tests  on  whole  aircraft  and  on  individual 
components.  Survivability/vulnerability  test  technology  has 
received  a  good  deal  of  attention  in  areas  such  as  test  specimen 
design,  threat  simulation,  and  simulation  of  the  aircraft  flow 
environment  seen  by  combat  aircraft.  One  aspect  cf  survivability/ 
vulnerability  testing  which  has  not  been  adequately  treated  in 
the  past,  is  the  simulation  of  actual  flight  loadings  of 
aircraft  structures  during  ballistic  impact  when  the  airflow 
environment  is  also  being  simulated. 

Another  way  of  assessing  the  survivability/vulnerability 
characteristics  of  aircraft  is  to  perform  large  scale  structural 
analyses.  In  the  past,  most  structural  analyses  of  aircraft 
have  been  restricted  to  linear,  small  deformation  applications. 
Such  analyses  have  not  shown  particularly  good  agreement  with 
experimentally  obtained  data  when  applied  to  ballistically 
damaged  structures,  however.  In  order  to  predict  load 
redistributions  and  residual  strength  in  damaged  structures, 
nonlinear  analyses  which  account  for  elastic-plastic  effects 
and  large  deformations  may  be  necessary.  These  extended  range 


analyses  have  been  applied,  to  a  limited  extent,  to  damaged 
aircraft  structures,  but  they  are  usually  very  inconvenient  for 
a  survivability/vulnerability  engineer  to  use  effectively. 

This  report  describes  the  results  of  a  research  project 
which  was  intended  to  investigate  the  possibility  of  removing 
some  of  the  shortcomings  in  both  the  experimental  and  analytical 
phases  of  survivability/vulnerability  studies.  The  objectives 
of  the  study  are  indicated  below. 

1.1  OBJECTIVES 

The  primary  objective  of  this  effort  was  to  develop  a 
validated  experimental  test  capability  for  the  realistic 
simulation  of  flight  loads  on  aircraft  wing  structures  during 
ballistic  impact.  In  achieving  this  primary  objective,  two 
equally  important  secondary  objectives  were  identified: 

a.  the  development  of  a  self-contained  flight  loads 
simulation  fixture, 

b.  the  development  of  a  structural  analysis  technique 
for  the  analytical  prediction  of  internal  load  distributions 
of  damaged,  multiple  load  path  structures. 

1.2  GENERAL  APPROACH 

The  general  approach  to  completing  the  program  objectives 
is  indicated  in  the  schematic  of  Figure  1.1.  Briefly,  the  total 
effort  is  split  initially  into  parallel  experimental  and 
analytical  developments.  These  experimental  and  analytical 
efforts  coalesce  into  an  error  analysis  of  data  obtained  from 
experimental  and  analytical  tests  using  common  replica  test 
specimens.  After  the  error  analysis  shows  comparable  results 
(perhaps  after  several  iterations  as  shown  by  the  double-headed 
arrows) ,  the  full-scale  test  item  is  mounted  in  the  experimental 
and  analytical  fixtures.  Then  reasonable  agreement  between  the 
experimental  and  analytical  test  results  for  the  full-scale 
item  is  taken  to  constitute  a  validation  of  the  experimental 
test  facility. 
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SECTION  2 

MAJOR  DEVELOPMENT  ITEMS 


This  section  describes  the  three  major  development  items: 
the  experimental  facility,  the  structural  analysis  technique, 
and  the  replica  test  specimens.  These  items  are  discussed  in 
more  detail  in  References  1-3;  the  brief  descriptions  presented 
here  are  included  for  completeness. 

2.1  EXPERIMENTAL  FACILITY 

The  experimental  facility’*'  designed  and  fabricated  by  the 
University  of  Dayton  Research  Institute  consists  of  a  self- 
contained  loading  fixture  and  a  loading  control  system.  These 
systems  impose  realistic  flight  loads  to  aircraft  wing  sections. 
When  mounted  in  the  Vertical  Gunfire  Facility  at  Wright-Patterson 
Air  Force  Base,  the  test  facility  applies  the  flight  loads 
simultaneously  with  ballistic  impact  and  airflow.  The  following 
paragraphs  present  discussions  of  the  basic  concept  of  the 
experimental  facility,  the  test  fixture,  and  the  control  system. 

2.1.1  Concept 

In  order  to  conduct  tests  with  realistic  flight 
loads,  fixtures  have  been  required  in  the  past  which  have  been 
large  enough  to  mount  whole  aircraft  wings;  or  tests  on  complete 
aircraft  have  been  performed.  The  objective  in  this  project 
was  to  design  an  experimental  facility  which  would  apply 
realistic  flight  loads  to  wing  sections  small  enough  so  that 
the  tests  could  be  performed  inside  the  Vertical  Gunfire  Range 
at  Wright-Patterson  Air  Force  Base. 

Figure  2.1  shows  diagrammatically  a  complete 

wing  and  a  section  cut  from  the  wing.  The  isolated  wing  section 

has  certain  values  of  spanwise  bending  moment  (M  ) ,  spanwise 

s 

shearing  force  (Vg)  chordwise  bending  moment  (Mc) ,  chordwise 
shearing  force  (V  ) ,  and  torque  (T) ,  acting  on  the  ends  of  the 


section.  The  values  of  the  section  loads  depend  on  the  loads 
applied  to  the  wing  in  specific  instances.  The  function  of  the 
experimental  facility  is  to  subject  wing  sections  to  operator 
specified  values  of  the  section  loads,  Mg ,  Vgf  Mc,  Vc#  and  T. 

Figure  2.2  shows  conceptually  the  experimental 
facility  which  has  been  designed,  fabricated,  and  installed  at 
Wright-Patterson  Air  Force  Base.  The  facility  generally 
consists  of: 

•  a  test  fixture  to  apply  loads  to  one  end  of 
a  specimen  and  to  react  at  the  other  end; 

•  a  hydraulic  system  to  impose  on  the  loading 
frame,  actuator  forces  which  correspond  to  specified  values 
of  the  section  loads;  and 

•  a  control  system  to  provide  the  test  operator 
with  a  convenient  means  for  controlling  the  application  of  the 
actuator  loads  during  a  test. 

A  test  on  a  wing  section  is  performed  by  first 

mounting  one  end  of  the  specimen  to  the  reaction  fixture  which 

is  designed  to  be  sufficiently  massive  and  rigid  so  that 

negligible  support  movement  occurs.  Loads  are  then  applied  to 

the  other  end  of  the  specimen  through  a  rigid  loading  frame. 

Strategically  placed  and  oriented  hydraulic  actuators  impose 

forces  on  the  loading  frame  which  in  turn  imposes  section 

loads  (M  ,  V  ,  M  ,  V  ,  T)  to  the  specimen.  The  loads  in  the 
s  s  c  c 

actuators  are  controlled  by  varying  the  pressure  in  the 
hydraulic  lines  through  a  remote  control  console.  The  control 
console  provides  the  capability  for  adjusting  the  relative 
magnitudes  of  the  actuator  loads  and  then  to  increase  them 
proportionally.  Note  that  the  flight  loads  can  be  applied 
simultaneously  with  air  flow  (existing  in  the  Air  Force 
facility)  and  ballistic  impact  (existing  in  the  Air  Force 
facility) . 

Figure  2.3  shows  a  damaged  wing  section  mounted 
in  the  experimental  facility.  A  typical  test  on  a  damaged 
specimen  consists  of  proportionally  increasing  the  actuator 
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Figure  2.2.  Test  of  Undamaged  Wing  Section. 


Figure  2.3.  Test  of  Damaged  Wing  Specimen. 


loads  until  failure,  in  order  to  obtain  an  indication  of  the 
residual  strength.  If  the  proportional  loads  correspond  to  a  mix 
of  bending  moments,  shears,  and  torque  for  a  certain  flight 
condition,  then  the  failure  load  can  be  related,  for  example,  to 
a  certain  maximum  g-condition  for  the  maneuver. 

2.1.2  Test  Fixture 

The  test  fixture  (Reference  1)  designed  to  apply 
realistic  flight  loads  to  aircraft  wing  sections  is  a  self-contained 
structure.  The  only  load  imposed  on  the  structure  of  the  Vertical 
Gunfire  Facility  during  a  test  is  the  dead  weight  of  the  test 
fixture,  specimen,  and  instrumentation.  That  is,  the  test  fixture 
is  free-standing;  it  does  not  have  to  be  installed  in  the  Vertical 
Gunfire  Facility  if  airloads  and/or  ballistic  impact  are  not 
considered  during  a  test. 

Figures  2.4a-c  are  photographs  of  the  test  fixture 
taken  from  various  angles.  The  photographs  show  all  of  the  major 
parts  of  the  fixture: 

•  the  reaction  structure, 

•  the  loading  frame, 

•  the  base  structure, 

•  the  hydraulic  actuators. 

These  systems  together  with  a  test  specimen  form  a 
closed  structural  loop. 

The  reaction  structure  consists  of  a  thick  steel 

plate  mounted  on  a  massive  back-up  structure  to  eliminate  rigid 

body  motion  of  the  test  specimen.  One  end  of  a  test  specimen  is 

attached  to  the  reaction  structure  through  custom  designed  mounting 

brackets?  The  reaction  structure  is  designed  to  react  maximum 

loads  of  M  =  12,000,000  in.  lb.,  V  =  100,000  lb.,  M  =  1,440,000 
s  sc 

in.  lb.,  Vc  =  12,000  lb.,  and  T  =  2,000,000  in.  lb.  The  reaction 
structure  is  fabricated  from  corrosion  resistant,  high-strength,  low 
alloy  (ASTMA  588)  steel,  having  a  strength  of  50,000  psi. 

*The  mounting  brackets  proved  to  be  a  weak  link  in  anchoring 
the  test  specimens  to  the  back-up  structure.  See  the  discussion 
in  Section  3.2.1,  p  32. 


The  loading  frame  is  a  rigid,  built-up  plate 
which  transfers  concentrated  forces  from  the  hydraulic 
actuators  into  the  test  specimen.  The  specimen  is  attached 
to  the  loading  frame  through  custom  designed  mounting  brackets. 

The  base  structure  is  a  framework  of  I-beams  which 
supports  the  reaction  structure  at  one  end  and  the  hydraulic 
actuators  at  the  other  end.  The  support  for  the  actuators  is 
moveable  to  allow  for  different  length  specimens.  The  base 
structure  is  fabricated  from  corrosion  resistant  high-strength, 
low  alloy  steel. 

The  hydraulic  actuators  are  reacted  by  the  base 
structure  at  one  end  and  impose  forces  on  tne  loading  frame 
at  the  other  end.  The  general  arrangement  of  the  actuators 
is  shown  in  Figure  2.5.  Four  actuators  are  oriented  vertically 
and  one  actuator  is  oriented  horizontally  in  a  chordwise 
direction . 

2.1.3  Control  System 

Control  of  the  hydraulic  loading  system  is 
accomplished  by  electrical  command  signals  from  an  operator's 
control  console.  Hydraulic  pressure  regulators,  driven  by  a 
controlled  electric  current,  regulate  the  oil  pressure  supplied 
to  the  loading  actuators.  Load  cells  sense  the  resulting 
loads  and  give  the  console  operator  a  positive  readout  of  the 
loads  being  applied.  The  operator's  console  contains  provisions 
for  setting  the  ratios  between  the  various  actuator  loads  and 
also  for  controlling  all  loads  in  unison.  A  block  diagram 
of  the  control  scheme  is  shown  in  Figure  2.6. 

The  control  console  (Figure  2.7)  contains  all 
controls,  readouts,  and  other  interfaces  with  the  operator, 
and  all  electronic  equipment  associated  with  the  loading 
system.  The  console  is  a  free  standing,  50  in.  high,  24  in. 
wide  cabinet  with  sloping  control  panel  (Figure  2.8).  Front 
panel  controls  include  the  key-operated  master  switch,  hydraulic 
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- *  ELECTRICAL  CONNECTIONS 
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- PILOT  RETURN  OIL 


Figure  2.6.  Block  Diagram  of  Load  Control  System 


LOAD  SYSTEM  CONTROL  AND  MONITORING 
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power  unit  control,  emergency  dump,  master  load  control,  and 
five  individual  load  controls.  Deflection  limit  set  dump 
controls,  calibration  span  set  controls,  and  calibration 
switches  also  are  included.  Digital  readouts  include  master 
load  control  setting,  five  load  readouts,  system  hydraulic 
pressure,  and  three  specimen  deflections. 

2.2  STRUCTURAL  ANALYSIS  TECHNIQUE 

2 

The  structural  analysis  technique  provides  the 
survivability/vulnerability  engineer  with  an  analytical  tool 
for  predicting  the  static  response  of  undamaged  and  damaged 
wing  structures.  The  objectives  of  the  structural  analysis 
capability  are:  (a)  to  provide  a  check  on  experimentally 
obtained  data,  (b)  to  compute  the  internal  stresses  in  undamaged 
wing  structures,  (c)  to  predict  the  stress  redistribution  in 
damaged  wing  structures,  and  (d)  to  estimate  the  residual 
strength  of  damaged  wing  structures. 

The  general  form  of  the  structural  analysis  technique 
is  indicated  in  Figure  2.9.  This  total  analysis  package 
consists  of  three  major  blocks  of  computer  code.  The  primary 
part  of  the  analysis  tool  is  a  finite  element  structural 
analysis  computer  program,  while  the  other  two  programs  are  a 
preprocessor  and  a  postprocessor.  The  pre-  and  postprocessors 
essentially  "straddle"  the  finite  element  program  in  a  manner 
so  that  the  composite  analysis  tool  is  convenient  for  the 
survivability/vulnerability  engineer  to  use  effectively. 

The  programs  shown  in  Figure  2.9  are  discussed  briefly 
in  the  following  paragraphs.  Detailed  instructions  for  using 
the  structural  analysis  technique  are  contained  in  Reference  2. 

2.2.1  Finite  Element  Program 

The  finite  element  program  is  the  major  part  of  the 
structural  analysis  technique.  The  particular  program  used  is 
the  MAGNA  (Materially  And  Geometrically  Nonlinear  Analysis) 


18 


PREPROCESSOR 


Figure  2.9 


Structural  Analysis  Technique. 
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finite  element  program  .  This  program  was  initially  developed 
by  the  University  of  Dayton  Research  Institute  using  internal 
funds.  Therefore,  MAGNA  is  considered  proprietary,  although 
only  to  a  limited  extent  as  far  as  the  Air  Force  is  concerned. 
Portions  of  the  later  stages  of  development  of  MAGNA  have  been 
funded  by  Air  Force  contracts.  Therefore,  the  only  restriction 
placed  on  the  MAGNA  program  for  the  Air  Force  is  that  it  not 
be  distributed  to  any  non-Air  Force  organization  without 
permission  from  the  University.  MAGNA  can  be  utilized  on  the 
Wright-Patterson  ASD  CDC  computers  by  both  Air  Force  personnel 
and  by  Air  Force  contractors.  Refer  to  proprietary  data 
statement  included  in  Reference  4. 

MAGNA  is  a  large-scale  computer  program  for  the  static 
and  dynamic  analysis  of  complex,  three-dimensional  engineering 
structures.  The  program  is  based  upon  the  finite  element  method 
of  analysis  to  permit  the  simulation  of  practical  structures 
composed  of  many  different  types  of  elements.  MAGNA  combines 
effective  isoparametric  modeling  techniques  with  state-of-the- 
art  numerical  analysis  and  programming  methods  to  provide 
accurate  and  efficient  solutions  for  large  problems  involving 
highly  nonlinear  response. 

The  modeling  capabilities  of  MAGNA  include  structural 
elements  for  truss  members,  plane  stress  and  plane  strain 
sections,  "shear  panels,"  general  three-dimensional  solids,  and 
thin  plates  and  shells.  All  finite  elements  are  arbitrarily 
oriented  and  are  fully  compatible  in  three-dimensional  space. 
Degrees  of  freedom  can  be  coupled  to  represent  skewed 
boundary  conditions,  rigid  regions,  and  complex  structural 
joints.  Uniform  mass  damping,  as  well  as  structural  damping 
based  upon  the  instantaneous  stiffness,  can  be  applied  in  the 
solution.  Time  history  solutions  are  performed  in  MAGNA 
using  an  implicit  scheme  for  direct  integration  of  the 
equations  of  motion. 
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Each  of  the  finite  elements  in  MAGNA  includes  the  effects 


of  full  geometrical  nonlinearities  (large  displacements, 
large  strains) ,  using  a  Lagrangian  (fixed  reference)  description 
of  motion.  In  shell  analysis,  arbitrarily  large  rotations  can 
also  be  treated.  Material  nonlinearities,  in  the  form  of 
elastic-plastic  behavior,  are  analyzed  using  a  subincremental 
strategy  which  minimizes  the  error  in  following  the  material 
stress-strain  curve.  Isotropic,  kinematic,  and  combined 
strain-hardening  rules  are  available  for  use  in  plastic  analysis 
with  MAGNA. 

The  MAGNA  program  includes  numerous  user  convenience 
features  to  aid  in  the  generation  of  finite  element  modeling 
data.  Geometry  data  may  be  input  in  Cartesian,  cylindrical, 
and  spherical  coordinates,  or  in  arbitrary,  user-defined 
systems.  Incremental  generation  of  nodal  coordinates  and 
element  connections  is  also  available  to  exploit  repetitive 
patterns  in  the  structural  model.  User-written  subroutines, 
which  provide  for  user  intervention  or  specif ication  of  data 
at  several  stages  of  the  analysis,  can  be  supplied  for 
defining  mesh  geometry,  coordinate  systems,  and  incremental 
applied  loading. 

Plotting  utilities,  in  both  interactive  and  batch  forms, 
are  also  available  for  use  in  checking  data,  and  for  interpreting 
analysis  results  obtained  from  MAGNA.  Geometry  plotting, 
including  exploded  views,  is  currently  available  for  all  finite 
elements.  Postprocessing  functions,  which  are  presently 
provided  for  most  element  types,  include  stress  and  strain 
contours  and  stress  relief  plots.  Scaled  and  exploded  views 
of  close-up  plots  of  the  deformed  structural  model  can  be 
generated,  with  the  undeformed  geometry  optionally  superimposed 
in  the  display. 

2.2.2  Preprocessor  Program 

The  preprocessor  program,  WINGEN,  is  a  convenient 
and  flexible  computer  program  designed  specifically  for  the 
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generation  of  wing  finite  element  models.  The  WINGEN  program 
accepts  data  in  an  abbreviated  format  and  then  generates  an 
expanded  set  of  data  which  is  acceptable  by  the  finite  element 
program. 

The  purpose  of  the  pre-processor  is  to  provide 
the  survivability/vulnerability  engineer  with  the  capability 
to  use  a  powerful,  state-of-the-art,  nonlinear  finite  element 
analysis  computer  program  without  having  to  spend  the  time 
and  effort  necessary  for  manual  preparation  of  detailed  wing 
finite  element  models  and  the  associated  data  decks.  Typical 
data  which  is  input  to  WINGEN  includes  coordinates  of  key 
points  to  define  the  planform  outline  of  the  wing,  number  and 
spacing  of  spars,  number  and  spacing  of  ribs,  skin  and  rib 
thicknesses,  bar  areas  and  material  properties.  The 
pre-processor  assumes  that  the  input  data  refers  to  a 
predetermined  class  of  structures  (in  particular,  wing 
structures).  The  preprocessor  automatically  completes  a 
finite  element  model  of  the  wing,  including  the  numbering  of 
the  nodes  and  the  elements.  Then  data  for  the  MAGNA  finite 
element  program  is  generated  in  the  appropriate  format;  typical 
data  generated  are  node  coordinates,  element  connectivities, 
degree  of  freedom  numbers,  properties  data,  and  program  control 
parameters. 

The  preprocessor  assumes  that  the  wing  skins 
are  two-dimensional  membrane  finite  elements,  the  spar  and 
rib  webs  are  two-dimensional  shear  panel  finite  elements,  and 
the  spar  and  rib  caps  are  one-dimensional  bar  finite  elements. 

That  is,  individual  skin  and  web  panels  do  not  resist  local 
bending  forces;  therefore,  the  resulting  finite  element  model 
for  some  damage  cases  may  be  unrealistic.  This  restriction 
was  defined,  however,  in  a  meeting  with  the  Air  Force.  One 
area  for  future  development  on  the  wing  preprocessor  is  to 
expand  the  types  of  finite  elements  which  can  be  automatically 
generated  in  wing  models.  Note  that  the  restriction  of  finite 
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element  types  to  bars,  membranes,  and  shear  panels  refers 
only  to  the  preprocessor;  the  MAGNA  program  can  solve  more 
sophisticated  models,  but  at  present  the  data  would  have  to 
be  prepared  manually. 

The  preprocessor  also  has  the  capabilitv 
to  create  finite  element  models  of  wing  structures  which  are 
damaged.  The  WINGEN  preprocessor  creates  damage  models  in 
two  ways : 

1.  The  user  defines  a  list  of  element  numbers, 
and  those  elements  are  dropped  from  the  undamaged  model,  and 
a  new  model  is  generated;  and 

2.  The  user  defines  the  coordinates  of  the 
center  of  a  sphere  and  a  radius,  and  the  preprocessor  drops 
from  the  undamaged  model  all  elements  having  centroids  within 
the  sphere,  and  generates  a  new  model. 

The  output  of  the  preprocessor  is  a  file  containing 
the  data  for  the  finite  element  program.  This  output  file  is 
also  readable  by  the  postprocessor  so  that  undeformed  geometry 
plots  can  be  viewed  before  running  the  finite  element  program. 

2.2.3  Postprocessor 

The  postprocessing  capability  consists  of  two 
computer  programs; 

1.  PLOTBOB  -  This  interactive  postprocessor 
accepts  data  files  from  either  the  preprocessor  or  from  the 
finite  element  program.  Data  from  the  preprocessor  is  used  to 
present  plots  of  undeformed  finite  element  meshes  for  data 
checking  purposes.  Data  from  the  finite  element  program  is 
used  to  present  plots  of  the  deformed  structure  either  singly 
or  superimposed  on  top  of  the  undeformed  structure.  Some  of 
the  options  of  this  postprocessor  are;  plotting  of  lists  of 
elements,  zooming,  clipping,  and  exploding.  The  model  can 
be  rotated,  translated,  or  reflected  about  defined  axes.  Also, 
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either  orthographic  projection  or  perspective  views  can  be 
selected.  In  addition,  the  user  can  selectively  label  nodes 
and  elements. 

2.  CONTOUR  -  This  interactive  postprocessor  is 
designed  to  display  the  results  of  the  finite  element  program 
in  various  forms  which  the  user  requests  selectively.  Examples 
of  the  plotting  options  available  are: 

•  contour  plots  of  stress,  strain,  and 
displacement  on  either  the  undeformed  or  the  deformed  model, 

•  selective  labeling  of  contours, 

•  relief  maps  of  stress,  strain  and  displacement, 

•  undeformed  and  deformed  geometry, 

•  zooming,  clipping,  and  exploding, 

•  rotation,  translation,  and  reflection,  and 

•  descriptive  plot  labeling. 

2.3  REPLICA  SPECIMENS 

The  replica  specimens  were  designed  to  demonstrate  and 
evaluate  the  performance  of  the  experimental  facility  and 
the  structural  analysis  technique.  The  specimens3  are 
simple,  rectangular  construction  to  allow  comparison  of 
experimental  and  analytical  results  without  introducing 
possible  modeling  errors  due  to  geometrical  complexity. 

A  photograph  of  a  typical  replica  test  specimen  is  shown 
in  Figure  2.10.  The  specimens  are  aluminum  having  overall 
dimensions  of  18  in.  deep,  60  in.  wide,  and  94.5  in.  long. 

The  cross-section  of  the  four-spar  specimens  is  shown  in 
Figure  2.11,  and  the  planform  is  shown  in  Figure  2.12. 

The  total  set  of  six  replica  specimens  consists  of  two 
undamaged  and  four  damaged  articles.  Figure  2.12  indicates 
a  scheme  for  referring  to  individual  skin  panels  (e.g.,  2-4), 
spar  web  panels  (B-5) ,  rib  web  panels  (d-3),  spar  caps  (C-3), 
and  rib  caps  (e-2).  The  two  end  bays  are  covered  by  the 
mounting  brackets  (see  Figure  2.4c)  and  are  not  included  in 


24 


iquro  2.10.  Replica  Test  Specimen 


of  Replica  Test  Spec 


the  analysis  model.  Therefore,  for  convenience  of  the  analysis, 
the  end  bays  are  not  numbered. 

The  physical  characteristics  of  the  six  replica  specimens 
are  listed  in  Table  2.1.  Specimen  Number  1  had  no  damage  and 
was  used  as  a  control  specimen  to  perform  initial  checkout 
of  the  experimental  facility.  Specimens  2,  3,  4,  and  5  contained 
simulated  damage  as  indicated  in  Table  2.1.  Specimen  6  was  a 
backup  specimen.  The  damaged  specimens  were  intended  to  check 
out  the  the  capability  of  the  experimental  facility  to  apply 
realistic  flight  loads  to  damaged  wing  sections  and  produce 
failure  in  the  test  articles.  The  replica  specimens  also 
provided  experimental  results  to  compare  with  analytically 
determined  stresses  at  various  points  of  the  damaged  specimens. 


TABLE  2.1 
REPLICA  SPECIMENS 


Specimen 

Number  Damage 


1 

None 

2 

Skin 

,  2-4 

,  Bottom,  Removed 

3 

Skin 

1 

CM 

,  Bottom  and  Top,  Removed 

Spar 

Web, 

B-4, 

C-4 ,  Removed 

Spar 

Cap, 

B-4, 

C-4 ,  Bottom  and  Top, 

Removed 

4 

Skin 

,  1-3 

,  Bottom  and  Top,  Removed 

Spar 

Web, 

A-3 , 

Removed 

Spar 

Cap, 

A-3  , 

Bottom  and  Top,  Removed 

Spar 

Cap, 

B-4, 

B-5 ,  Bottom  and  Top, 

Split 

5 

Skin 

,  1-4 

,  2-4 

,  3-4,  Bottom,  Removed 

Spar 

Web, 

B-4  , 

C-4 ,  Removed 

Spar 

Cap, 

B-4, 

C-4,  Bottom  and  Top, 

Removed 

6 

None 
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SECTION  3 

REPLICA  SPECIMEN  TESTS  AND  ANALYSES 


This  section  discusses  the  results  of  tests  performed 
on  the  replica  specimens  using  the  experimental  facility,  and 
of  analyses  of  finite  element  models  of  the  replica  specimens 
using  the  structural  analysis  technique. 


3.1  ACTUATOR  LOADS  AND  EQUIVALENT  SECTION  LOADS 


In  the  tests  and  analyses,  the  section  loads  (spanwise 
bending  moment  Mg,  spanwise  shear  force  V  ,  torque  T,  and 
chordwise  shear  force  V  )  at  the  unsupported  end  of  a  test 
specimen,  are  assumed  to  be  known.  The  positions  of  the 
hydraulic  actuators  relative  to  the  loaded  end  of  the  replica 
specimens  are  shown  schematically  in  Figure  2.5.  Positive 
values  of  the  section  loads  and  the  actuator  forces  (tension 
positive)  are  shown  in  Figure  3.1.  Expressions  which  relate 
the  actuator  forces  and  the  applied  section  loads  are: 


T1  + 


t2  + 


T3  + 


Tlg  ~  T2^  ”  T3c  +  T4^  =  _T 

(3.1) 

T, b  +  T_b  -  T0a  -  T.a  =  -M 
1  2  3  4  s 

T4a  -  T]b  =  T3a  -  T2b 


T 


5 


-V 

c 


The  fourth  equation  ensures  that  the  actuator  forces  all  work 
approximately  equally;  the  equation  forces  the  bending  moment 
due  to  pairs  of  actuators  on  one  side  of  the  center  line  of 
a  specimen  to  be  equal  to  the  bending  moment  due  to  pairs  of 
actuators  on  the  other  side  of  the  center  line. 

Solving  Equations  3.1  for  the  actuator  forces  qives: 
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Figure  3.1.  Positive  Sign  Convention. 


T_  =  -V 
5  c 

These  equations  are  used  to  determine  what  actuator  loads  the 
operator  of  the  control  console  must  impose  in  order  to  apply 
given  values  of  the  section  loads. 

3.2  TEST  1  -  UNDAMAGED  SPECIMEN,  NUMBER  1 

The  first  test  was  a  series  of  verification  tests  to 
determine  if  the  experimental  facility  was  operating  properly. 

A  number  of  simple  load  cases  such  as  pure  spanwise  bending, 
pure  spanwise  shear,  pure  torque,  and  combinations  thereof 
were  performed  with  undamaged  specimen  Number  1  (Table  2.1). 

3.2.1  Instrumentation 

Figure  3.2a  shows  a  schematic  representation  of 
the  relative  locations  of  strain  gages  monitored  during  the 
tests.  The  numbers  in  circles  identify  strain  gage  rosettes 
while  the  uncircled  numbers  indicate  individual  gages. 

Specimen  Number  1  also  was  instrumented  for  measuring 
deflections  at  several  points  on  the  bottom  surface  of  the 
specimen.  No  useful  deflection  data  was  obtained,  however,  due 
to  excessive  deformation  in  the  support  brackets  which  were  used 
to  mount  the  test  specimen  to  the  reaction  fixture.  The  support 
bracket  deformation  resulted  in  erroneous  deflection  measurements. 
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In  addition,  the  deflection  measurements  were  not  repeatable. 
Consequently,  the  measurement  of  deflections  was  abandoned  for 
the  remainder  of  the  replica  specimen  tests. 

3.2.2  Load  Conditions 

The  particular  loading  conditions  applied  in  the 
initial  testing  phase  are  given  in  Table  3.1.  Using  Equations  3 
the  individual  actuator  loads  used  for  each  of  the  above  loading 
cases  were  calculated,  these  actuator  loads  are  shown  in  Table  3 

3.2.3  Load  Incrementation  and  Data  Collection 

During  the  five  verification  tests,  the  actuators 
were  increased  in  increments  of  20  percent  of  the  maximum  loads 
listed  above  until  the  maximum  values  were  attained;  then  the 
actuator  loads  were  decreased  in  20  percent  increments.  Strain 
gage  readings  were  recorded  at  every  increment  of  loading  and 
unloading  by  a  minicomputer.  A  small  computer  program  converted 
the  gage  signals  into  units  of  strain.  In  addition,  the 
individual  gage  stresses  were  calculated,  and  in  the  case  of 
the  rosettes  the  minimum  principal  stress,  the  maximum  principal 
stress,  the  maximum  shear  stress,  and  the  principal  angle. 

These  calculated  values  were  printed  for  each  loading  case  for 
each  loading  increment.  An  example  of  the  output  of  the 
program  is  shown  in  Figure  3.2b.  The  output  corresponds  to  40 
percent  of  the  maximum  loading  of  case  1,  pure  bending. 

3.2.4  Analysis  Model 

The  undamaged  replica  specimen  finite  element 
model  is  shown  in  Figure  3.3a-f.  The  portion  of  the 
replica  specimens  modeled  is  that  part  of  Figure  2.12  which 
have  the  bays  numbered.  As  mentioned  before,  in  the 
experimental  facility  the  two  end  bays  are  clamped  by  mounting 
brackets  and  are,  therefore,  not  considered  in  the  analysis. 

The  model  contains  56  nodes,  36  skin  membrane  elements,  24  spar 
web  shear  panel  elements,  15  rib  web  shear  panel  elements,  48 
spar  cap  bar  elements,  and  30  rib  cap  bar  elements.  Each  of  the 
56  nodes  has  three-degrees  of  freedom  -  the  displacements 
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TABLE  3.2 
ACTUATOR  LOADS 


Verification 

Case 

Actuator  Loads  (lb)  j 

T1 

T 

2 

T3 

1 

-21,252 

-23,010 

+18,889 

+25,380 

2 

-  1,704 

+  1,704 

+  6,294 

-  6,294 

3 

-  3,069 

-  3,324 

-12,273 

-11,334 

4 

-15,700 

-16,998 

+  6,446 

+11,253 

5 

-17,404 

-15,294 

+12,740 

+  4,959 
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Figure  3.2b.  Sample  Experimental  Data  -  Specimen 
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Figure  3.2b.  (concluded). 
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Figure  3.3a.  Specimen  1  Finite  Element  Model  -  Node  Numbers. 


Figure  3.3c.  Specimen  1  Finite  Element  Model  -  Upper  Skin  Membrane  Elements. 


Figure  3.3f.  Specimen  1  Finite  Element  Model  -  Cap  Bar  Elements. 


parallel  to  the  three  coordinate  axes.  Eight  of  the  nodes  are 
clamped  at  the  reaction  end  of  the  specimen.  Therefore,  the 
analysis  model  has  3  x  48  =  144  degrees  of  freedom. 

3.2.5  Test/Analysis  Results 

The  series  of  load  cases  defined  in  Paragraph 
3.2.2  was  used  to  determine  if  the  experimental  facility 
operated  as  it  was  designed.  The  maximum  load  levels  were 
intentionally  kept  relatively  low  to  ensure  that  no  failure 
of  the  specimen  would  occur  during  the  tests. 

Each  of  the  five  loading  cases  was  applied  to 
the  specimen  by  incrementing  the  individual  actuator  forces 
through  the  control  console.  The  repeatability  of  the  tests 
was  checked  by  cycling  through  the  tests  twice.  In  each  case 
the  repeatability  was  excellent.  The  symmetric  load  cases 
produced  symmetric  strain  gage  readings.  Also,  the  strains 
on  the  top  and  bottom  surfaces  of  the  specimen  at  the  same 
planform  locations  were  of  the  same  magnitude  but  of  opposite 
signs. 

The  comparison  between  the  results  obtained  from 
the  verification  tests  and  those  produced  by  the  MAGNA  finite 
element  program  was  quite  good.  As  mentioned  above  in  Paragraph 
2.2.2,  the  finite  element  model  produced  by  the  wing  model 
preprocessor  is  a  one-element-per-bay  model.  This  type  of 
model  has  been  used  by  many  aircraft  stress  analysts  for 
analyzing  undamaged  wing  structures  with  good  success. 

Therefore,  it  was  expected  that  the  test/analysis  comparison 
for  the  undamaged  specimen  tests  would  be  good. 

As  an  example,  Figure  3.4  shows  the  comparison 
between  a  stress  measured  during  the  test  of  Case  1  (pure 
bending) ,  and  the  stress  predicted  by  the  finite  element  program. 
The  location  of  the  compared  stresses  corresponded  to  rosette 
number  3  in  Figure  3.2  and  element  29  in  Figure  3.3b.  Figure 
3.4  is  a  plot  of  the  maximum  principal  stress  at  the  indicated 
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location  versus  the  percentage  of  the  maximum  loading  applied 
to  the  specimen.  The  figure  shows  that  the  test  results  and 
the  finite  element  results  are  very  close.  Also  plotted  on 
the  figure  is  the  result  obtained  from  the  simple  beam  formula 
a  =  Mc/I.  The  close  agreement  between  the  three  results  (test, 
finite  elements,  and  beam  theory)  gives  insight  into  why  the 
bar/membrane/shear  panel  finite  element  modeling  approach  works 
so  well  for  undamaged  wings.  That  is,  a  wing  box  (undamaged) 
responds  sufficiently  like  a  classical  beam  so  that  a  simplified 
finite  element  model  provides  accurate  analytical  results. 
However,  this  is  not  the  case  for  some  damaged  wing  structures 
(see  Sections  3. 3-3. 6). 

The  results  of  the  remaining  verification  tests 
were  similar  to  those  of  Case  1  (pure  bending)  discussed 
above.  Therefore,  it  was  concluded  that  the  experimental 
facility  operated  as  designed.  Some  minor  problems  with  the 
electrical  and  hydraulic  systems  were  encountered,  but  these 
were  considered  to  be  more  of  a  debugging  nature  rather  than 
design  deficiencies. 

3.3  TEST  2  -  DAMAGED  SPECIMEN,  NUMBER  2 

The  second  test  was  performed  on  a  specimen  having  a 
small  amount  of  damage  -  one  panel  missing  from  the  lower  skin. 

A  combination  of  spanwise  bending  and  spanwise  shear  loads 
were  applied.  The  intention  was  to  increment  the  applied  loads 
until  failure  of  the  specimen  occurred.  The  following 
paragraphs  describe  Test  2  in  detail. 

3.3.1  Instrumentation 

Figure  3.5  shows  a  schematic  representation  of 
the  relative  locations  of  strain  gages  monitored  during 
Test  2.  The  numbers  in  circles  identify  strain  gage 
rosettes,  while  the  uncircled  numbers  indicate  individual 


3.3.2 


Loading 


The  maximum  load  applied  to  a  specimen  was  a 
combination  of  spanwise  bending  moment  and  spanwise  shear  load 
(Mg  =  5.4  x  10^  in.  lb.  and  Vg  =  30,000  lb).  The  individual 
actuator  loads  corresponding  to  these  section  loads  were: 

T^  =  45,573  lb.  Compression 

T 2  =  49,344  lb.  Compression 

T^  =  25,497  lb.  Tension 

T^  =  29,426  lb.  Tension 

3.3.3  Load  Incrementation  and  Data  Collection 

During  the  test  the  actuator  loads  were 
increased  in  5  percent  increments  until  failure  of  the  specimen 
occurred.  At  each  stage  of  the  loading,  strain  gage  readings 
were  recorded  by  a  minicomputer.  A  computer  program  converted 
the  gage  signals  into  units  of  strain,  and  computed  the 
corresponding  gage  stresses.  In  the  case  of  the  rosettes,  the 
maximum  principal  stress,  the  maximum  shear  stress,  and  the 
principal  angle  were  calculated.  An  example  of  the  output  of 
the  data  reduction  program  for  Test  2  is  shown  in  Figure  3.6. 
The  output  shown  corresponds  to  35  percent  of  the  maximum 
loads  defined  above. 

3.3.4  Analysis  Model 

The  finite  element  model  of  Specimen  Number  2  is 
shown  in  Figure  3.7a-e.  The  modeled  portion  of  the  replica 
specimen  is  that  part  of  Figure  2.12  which  has  the  bays 
numbered.  As  mentioned  before,  in  the  experimental  facility, 
the  two  end  bays  are  clamped  by  mounting  brackets,  and  are 
therefore  not  considered  in  the  analysis.  The  model  contains 
56  nodes,  35  skin  membrane  elements,  24  spar  web  shear  panel 
elements,  15  rib  web  shear  panel  elements,  48  spar  cap  bar 
elements,  and  30  rib  cap  bar  elements.  Each  of  the  56  nodes 
has  three  degrees  of  freedom  -  the  displacements  parallel  to 
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Figure  3,6.  Sample  Experimental  Data  -  Test 
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Figure  3.0.  (concluded). 
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Figure  3.7a.  Specimen  2  Finite  Element  Model  Node  Numbers. 


Ficmre  3.7c.  Specimen  2  Finite  Element  Model  -  Spa r  Shear  Panel  Elements. 


Figure  3.7e.  Specimen  2  Finite  Element  Model  -  Cap  Bar  Elements. 
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the  three  coordinate  axes.  Eight  of  the  nodes  are  clamped 
at  the  reaction  end  of  the  specimen.  Therefore,  the  analysis 
model  has  3  x  48  =  144  degrees  of  freedom. 

3.3.5  Test/Analysis  Results 

The  replica  specimen  number  2  had  minimal  damage, 
with  a  single  skin  panel  missing  on  the  lower  skin  and  no  other 
structural  damage  present.  This  small  amount  of  damage  was  not 
expected  to  cause  the  response  of  the  specimen  to  deviate 
much  from  that  which  would  be  predicted  with  a  mathematical 
model  based  on  the  use  of  membranes,  shear  panels,  and  bars 
for  the  skins,  webs,  and  caps,  respectively.  Comparison  of  the 
experimental  and  analytical  results  showed  the  expectation  to 
be  true.  As  an  example,  Figure  3.8  shows  the  experimentally 
and  analytically  obtained  maximum  principal  stress  for  rosette 
number  8  (Figure  3.5)  and  finite  element  number  19  (Figure  3.7b). 
It  can  be  seen  that  the  analytical  prediction  is  quite  close 
to  the  experimental  result. 

The  loads  applied  to  the  specimen  were  incremented 
to  85  percent  of  the  maximum  loads  indicated  in  Paragraph  3.3.2. 
The  test  was  terminated  at  this  stage  because  of  excessive 
local  deformation  of  the  specimen  at  the  points  of  attachment 
to  the  reaction  structure  and  to  the  loading  frame.  The 
response  in  the  middle  of  the  specimen  remained  quite  linear, 
however,  as  seen  in  Figure  3.8. 

3.4  TEST  3  -  DAMAGED  SPECIMEN,  NUMBER  3 

Specimen  Number  3  (Paragraph  2.3)  had  more  extensive 
damage  than  Specimen  Number  2.  Figure  3.9  shows  a  schematic 
representation  of  the  damage  which  extends  the  entire  depth 
of  the  specimen.  Missing  are  skin  panels  on  the  upper  and 
lower  surfaces,  the  center  two  spar  shear  webs,  and  the  upper 
and  lower  spar  caps.  A  combination  of  spanwise  bending 
moment  and  spanwise  shear  loads  were  applied  to  the  specimen. 

The  following  paragraphs  describe  Test  3  in  more  detail. 
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3.4.] 


Instrumentation 


Figure  3.9  indicates  the  relative  locations  of 
strain  gages  monitored  during  Test  3.  Both  rosettes  and 
individual  gages  are  numbered. 

3.4.2  Loading 

The  maximum  load  applied  to  the  specimen  was 
a  combination  of  spanwise  bending  moment,  Mg  =  5.4  x  106  in.  lb. 
and  spanwise  shear  load,  Vg  =  30,000  lb.  The  individual 
actuator  forces  corresponding  to  these  section  loads  were 
(from  Equations  3.2): 

Tj  =  45,573  lb.  Compression 

T£  =  49,344  lb.  Compression 

T^  =  25,497  lb.  Tension 

T^j  =  29,426  lb.  Tension 

3.4.3  Load  Incrementation  and  Data  Collection 

During  the  test  the  .individual  actuator  forces 
were  increased  incrementally  in  steps  equal  to  5  percent  of 
the  maximum  values  given  above.  At  each  stage  of  the  loading, 
strain  gage  readings  were  recorded  by  a  minicomputer.  A 
computer  program  coverted  the  gage  signals  into  units  of 
strain,  and  computed  the  associated  material  stresses.  In 
the  case  of  the  rosettes,  the  minimum  and  maximum  principal 
stresses,  the  maximum  shear  stress,  and  principal  angle,  and 
the  Von  Mises  equivalent  stress  were  computed.  An  example 
of  the  output  of  the  data  reduction  program  for  Test  3  is 
shown  in  Figure  3.10.  The  output  shown  corresponds  to  60 
percent  of  the  maximum  loads  defined  above. 

3.4.4  Analysis  Model 

The  finite  element  model  for  specimen  No.  3 
is  shown  in  Figure  3.11a-e.  The  modeled  portion  of  the 
replica  specimen  is  that  part  of  Fioure  2.12  which  has  the 
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Figure  3.ila.  Specimen  3  Finite  Element  Model  -  Node  Numbers 


.lib.  Specimen  3  Finite  Element  Model  -  Bottom  Skin  Membrane  Elements. 


REPLICA  UINC  /  MOM- LI  ME  PR  ANALYSIS  /  DAMAGE  IH  CHORD  IRV  2,  SPAN  SAY  4 


Figure  3. lid.  Specimen  3  Finite  Element  Model  -  Spar  and  Rib  Shear  Panel  Elements. 


Figure  3. lie.  Specimen  3  Finite  Element  Model  -  Cap  Bar  Elements. 


bays  numbered.  As  mentioned  before,  in  the  experimental 
facility,  the  two  end  bays  are  clamped  by  mounting  brackets, 
and  are  therefore  not  considered  in  the  analysis.  The  model 
contains  56  nodes,  34  skin  membrane  elements,  22  spar  web  shear 
panel  elements,  15  rib  web  shear  panel  elements,  44  spar  cap  bar 
elemtns,  and  30  rib  cap  bar  elements.  Each  of  the  56  nodes  has 
three  degrees  of  freedom,  the  displacements  parallel  to  the 
three  coordinate  axes.  Eight  of  the  nodes  are  fixed  at  the 
reaction  end  of  the  specimen  (Nodes  1-8).  Therefore,  the 
analysis  model  has  3  x  48  =  144  degrees  of  freedom. 

3.4.5  Test/Analysis  Results 

When  the  console  operator  began  the  test  on 
Specimen  3,  a  malfunction  of  one  of  the  control  systems  caused 
one  of  the  hydraulic  actuators  to  overload.  The  result  of  the 
malfunction  and  overload  was  that  Specimen  3  was  damaged 
beyond  use  with  no  experimental  data  being  collected.  The 
spare  specimen  number  6  (Table  2.1)  was  then  altered  according 
to  the  damage  of  Specimen  3,  and  the  control  system  was 
modified  to  prevent  a  reoccurrence  of  the  malfunction.  The 
modified  specimen  6  is  referred  to  hereafter  as  Specimen  3. 

Figures  3.12a-p  compare  experimentally  obtained 
stresses  with  corresponding  stresses  computed  analytically  with 
the  finite  element  program.  Of  the  eight  skin  stresses  shown 
in  Figures  3.12a-h,  the  stresses  in  elements  19  and  21  (those 
elements  adjacent  to  the  damaged  area)  compared  most 
favorably;  Figures  3.13a-e  give  some  insight  into  why  this  is 
so.  Figures  3.13a-e  are  postprocessor  plots  which  present 
level  contours  of  the  equivalent  stresses  on  the  lower  skin 
superposed  on  the  finite  element  model  of  the  lower  skin.  In 
the  spanwise  bays  on  either  side  of  the  damaged  bay,  the 
stress  contours  undergo  considerable  change  as  the  stress 
"goes  around"  the  damaged  area.  Apparently,  the  "one-element- 
per-bay"  finite  element  model  is  not  sufficient  in  this  case 
to  predict  the  stress  redistribution  in t  he  transition  areas. 
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Figure  3.12d.  Equivalent  Stress-Skin  Element  19. 
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Figure  3.12h.  Equivalent  Stress-Skin  Element  27 


Figure  3.12k.  Equivalent  Stress-Spar  Cap  Element  35 
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Figure  3.12o.  Equivalent  Streas-Spar  Cap  Element  69 
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Figure  3.13e.  Contours  of  Equivalent  Stress 


The  discontinuity  which  appears  in  the 
experimental  data  in  Figures  3.12a-p  occurred  when  two  top 
skin  panels  (Elements  20  and  22  in  Figure  3.11c)  buckled. 

The  buckling  of  the  panels  caused  a  readjustment  in  the 
stresses  at  other  locations  on  the  specimen.  The  simplified 
finite  element  model  used  is  not  capable  of  predicting  local 
panel  buckling.  Also,  the  nonlinear  and  somewhat  erratic 
behavior  of  the  experimental  results  at  the  higher  load  levels 
is  attributed  not  only  to  the  compression  skins  buckling  but 
also  to  the  observed  failure  of  rivets.  The  simplified  finite 
element  model  also  does  not  predict  rivet  failure. 

3.5  TEST  4  -  DAMAGED  SPECIMEN,  NUMBER  4 

Specimen  4  is  unsymmetrical  having  a  "notch"  in  the 
leading  edge  and  longitudinally  split  interior  spar  caps  as 
indicated  in  Table  2.1.  Again  a  combination  of  spanwise 
bending  moment  and  spanwise  shear  load  were  applied  to  the 
specimen. 

3.5.1  Instrumentation 

Figure  3.14  shows  the  relative  locations  of 
strain  gages  monitored  during  Test  4.  Both  rosettes  and 
individual  gages  are  numbered  in  the  figure. 

3.5.2  Loading 

The  maximum  load  applied  to  the  specimen  was  a 
combination  of  spanwise  bending  moment  (Mg  =2.7  x  10^  in.  lb.) 
and  spanwise  shear  load  (V  =  30,000  lb.).  The  individual 
actuator  forces  corresponding  to  these  section  loads  were 
(from  Equations  3.2): 

T^  =  45,573  lb.  Compression 

T 2  =  49,344  lb.  Compression 

T^  =  25,497  lb.  Tension 


T.  =  29,426  lb.  Tension 


5/  f H'.tt  n  a 


Figure  3.14.  strain  Gage  Locations  -  Specimen 


3.5.3  Load  1  nc  remen  ta  t  ion  and  Data  Collection 

During  the  test  of  Specimen  4,  the  individual 
actuator  forces  were  increased  incrementally  in  steps  equal  to 
five  percent  of  the  maximum  values  qi ven  above.  At  each  stage 
of  the  loading,  strain  uaqe  readings  were  recorded  by  a  mini¬ 
computer.  A  computer  program  converted  the  gage  signals  into 
units  of  strain  and  computed  the  associated  material  stresses. 

In  the  case  of  the  rosettes,  t he  minimum  and  maximum  principal 
stresses,  the  maximum  shear  stress,  the  principal  angles,  and 
the  Von  Mises  equivalent  stress  were  computed.  An  example  of  the 
output  of  the  data  reduction  program  for  Test  4  is  shown  >n 
Figure  3.15.  The  output  shown  corresponds  to  40  percent  of  the 
maximum  loads  defined  above. 

3.5.4  Analysis  Model 

The  finite  element  model  for  specimen  4  is  shown  in 
Figure  3.16a-e.  The  modeled  portion  of  the  replica  specimen  is 
that  part  of  Figure  2.12  which  has  the  bays  numbered.  As 
mentioned  before,  in  the  experimental  facility,  the  two  end  bays 
of  the  specimen  are  clamped  by  mounting  brackets  and  are,  therefore, 
not  considered  in  the  analysis.  The  model  contains  56  nodes,  34 
skin  membrane  elements,  23  spar  web  shear  panel  elements,  15  rib 
web  shear  panel  elements,  46  spar  cap  bar  elements  and  30  rib 
cap  bar  elements.  Fach  of  the  56  nodes  has  three  degrees  of 
freedom,  the  disol  icemen t. s  parallel,  to  the  three  coordinate  axes. 
Fight  of  the  nodes  are  fixed  at  the  reaction  end  of  the  specimen 
(nodes  1-8).  Therefore,  the  analysis  model  has  3  x  48  -  144 
degrees  of  freedom. 

3.5.5  Test  7 Anal ysi s  Results 

Figures  i. I7a-m  compare  experimental ly  obtained 
stresses  with  correspond ing  stresses  computed  analytically  with 
the  finite  element  program.  Fach  plot  has  two  experimental 
results  displayed.  The  strain  qages  on  the  specimen  were  zeroed 
with  the  loadinq  apparatus  hanging  from  the  specimen.  The  lower 
curve  was  obtained  assuming  that  the  gages  had  zero  readings  when 
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Figure  3.14.  Sample  Experimental  Data  -  Specimen  4. 
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Figure  3.14.  (continued) 
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Figure  3.14  (concluded) 
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Figure  3.16d.  Specimen  4  Finite  Element  Model  -  Spar  and  Rib  Shear  Panel  Elements. 
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Figure  3.17m.  Equivalent  Stress-Spar  Web  Element  70. 
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the  actuator  forces  were  zero.  The  upper  curve  was  obtained  by 
takinq  into  account  the  initial  shear  and  moment  applied  to  the 
specimen  due  to  the  weiqht  of  the  loading  plate. 

There  was  rather  qood  comparison  between  the 
experimentally  obtained  and  analytically  computed  stresses  of 
skin  elements  13  and  15  (Figure  3.17a  and  b)  ,  the  elements 
adjacent  to  the  damaged  section  (Figure  3.16b).  There  was  only 
fair  agreement,  however,  for  the  other  skin  elements.  This 
again  indicates  that  the  bar/membrane/shear  panel  finite  element 
model  is  not  sufficient  to  predict  accurately  the  redistribution 
of  stresses  around  a  damaged  area.  Figures  3.17i-m  refer  to  shear 
webs.  The  analytical/experimental  comparisons  for  the  shear  web, 
run  from  fair  (Figures  3.17i,  1,  m,  n)  to  poor  (Figures  3.17j,  k) . 
The  erratic  behavior  of  the  gages  attached  to  elements  41  and  52 
indicate  that  there  might  have  been  some  malfunction  of  these 
channels,  however.  The  experimentally  measured  data  showed  a 
sudden  change  in  the  rate  of  stress  growth  at  about  the  70  percent 
load  level.  This  was  due  to  buckling  of  the  compression  skin 
around  the  area  of  the  split  spar  caps.  The  bar/membrane/shear 
panel  model  does  not  capture  this  buckling  effect. 

Figures  3.18a-e  present  plots  of  level  contours  of 
the  Von  Mises  equivalent  stress  at  various  stages  of  the  load 
incrementation . 

3.6  TEST  5  -  DAMAGED  SPECIMEN,  NUMBER  5 

Specimen  5  had  the  largest  amount  of  damage  with  all  of 
the  bottom  skin  and  all  of  the  interior  spar  webs  and  spar  caps  in 
one  bay  of  the  specimen  missino.  The  same  combination  of  spanwise 
bending  and  spanwise  shear  was  applied  to  this  specimen  as  was 
applied  to  the  other  damaged  specimens.  The  loading  was  increased 
proportionally  until  failure  of  the  specimen  occurred, 

3.6.1  Instrumentation 

Figure  3.19  shows  the  relative  locations  of  strain 
gages  monitored  during  Test  5.  Both  the  rosettes  and  individual 
gages  are  numbered  in  the  figure. 
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Figure  3.18b.  Contours  of  Equivalent  Stress  -  40%  Load 
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3.6.2 


Lo  ad  ing 


The  maximum  load  set  up  at  the  control  console 
corresponded  to  a  combination  of  snanwise  bendinq  moment 
(Ms  =  5.4  x  10^  in.  lb.)  and  spanwise  shear  load  (V  =  30,000  lb.). 
The  individual  actuator  loads  corresponding  to  these  section 
loads  were  (from  Equations  3.2): 

T^  =  45,573  lb.  Compression 

T2  -  49,344  lb.  Compression 

T^  =  25,497  lb.  Tension 

T^  =  29,426  lb.  Tension 

3.6.3  Load  Incrementation  and  Data  Collection 

During  the  test  of  Specimen  5,  the  individual 
actuator  forces  were  increased  incrementally  in  steps  equal  to 
5  percent  of  the  maximum  values  given  above.  At  each  stacie  of 
the  loading,  strain  gage  readings  were  recorded  by  a  minicomputer. 

A  computer  program  converted  the  gage  signals  into  units  of 
strain,  and  computed  the  associated  material  stresses.  In  the 
case  of  the  rosettes,  the  minimum  and  maximum  principal  stresses, 
the  maximum  shear  stress,  the  principal  angle,  and  the  Von 
Mises  equivalent  stress  were  computed.  An  example  of  the  output 
of  the  data  reduction  program  for  Test  5  is  shown  in  Figure  3.20. 
The  output  corresponds  to  20  percent  of  the  maximum  loads 
defined  above. 

3.6.4  Analysis  Model 

The  finite  element  model  for  specimen  5  is  shown 
in  Figure  3.21a-e.  The  modeled  portion  of  the  replica  specimen 
is  that  part  of  Figure  2.12  which  has  the  bays  numbered.  As 
mentioned  before,  in  the  experimental  facility  the  two  end  bays 
of  the  specimen  are  clamped  by  mounting  brackets,  and  are 
therefore  not  considered  in  the  analysis.  The  model  contains 
56  nodes,  33  skin  membrane  elements,  22  spar  web  panel  elements. 

15  rib  web  shear  panel  elements,  44  spar  cap  bar  elements,  and 
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Figure  3.21a.  Specimen  5  Finite  Element  Model  -  Node  Numbers 
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Specimen  5  Finite  Element  Model  -  Top  Skin  Membrane  Elements. 


REPLICA  UINC  /  NON-UNEAR  ANALYSIS  /  DAF1AGE  IN  CHORD  SAV  1-3  SPAN  IAV 


Specimen  5  Finite  Element  Model  -  Spar  and  Rib  Shear  Panel  Elements. 


30  rib  cap  bar  elements.  Each  of  the  56  nodes  has  three  degrees 
of  freedom,  the  displacements  parallel  to  the  three  coordinate 
axes.  Eight  of  the  nodes  are  fixed  at  the  reaction  end  of  the 
specimen  (nodes  1-8).  Therefore,  the  analysis  model  has 
3  x  48  =  144  degrees  of  freedom. 

3.6.5  Test/Analysis  Results 

Figures  3.22a-gg  compare  experimentally  determined 
stresses  with  corresponding  stresses  computed  with  the  finite 
element  program.  The  skin  element  stresses  are  considered  in 
Figure  3.22a-i.  The  best  comparison  between  experimental  and 
analytical  determinations  of  skin  stresses  are  for  elements 
19,  20,  21  (Figure  3.21c),  the  elements  on  the  undamaged  top 
surface.  For  skin  elements  13,  15,  17,  22,  24,  26,  the 
comparisons  are  not  so  close.  This  lack  of  comparison  is  again 
attributed  to  the  inability  of  the  simplified  bar/membrane/shear 
oanel  ^inite  element  model  to  capture  the  redistribution  of 
stresses  about  a  damaged  area.  Figures  3.23a--d  show  plots  of 
level  contours  of  equivalent  stress  in  the  bottom  skin  for 
various  load  levels.  All  of  the  stress  in  the  damaged  skin 
must  be  directed  around  the  damaged  area  through  the  spar  webs 
and  spar  caps  which  bound  the  damaged  area.  Apparently  the 
simplified  finite  element  model  is  too  coarse  to  predict  the 
stress  redistribution  accurately. 

The  stresses  in  rib  web  panels  40-45  are  considered 
in  Figures  3.22j-o.  Except  for  panel  43  the  analytical  and 
experimental  results  do  not  agree  well.  The  stresses  compare 
quite  well  in  the  case  of  the  spar  web  panels  51,  52,  53,  67,  66, 
69.  The  agreement  is  especially  good  for  loading  to  about  the 
70  percent  levels;  above  70  percent  the  two  results  deviated 
substantially  in  many  cases.  This  deviation  at  the  higher  load 
levels  is  attributed  to  the  occurrence  of  local  phenomena  such 
as  the  failure  of  rivets  and  the  buckling  of  compression  panels 
which  were  observed  during  the  test. 
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Figure  3.22b.  Equivalent  Stress-Skin  Element  15, 
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Figure  3.22c.  Equivalent  Stress-Skin  Element  17 
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Fiqure  i.22e.  Equivalent  Stress-Skin  Element  20 


T 


CO 

CO 

w 

cc 

Eh 

CO 


STRESS  VS  PERCENT  LORD-EL #2 4 


PERCENT  LOAD 


Figure  3.22h.  Equivalent  Stress-Skin  Element  24. 
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Figure  3.22t.  Equivalent  Stress-Spar  Web  Element  68. 
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Figure  3.22v.  Axial  Stress-Spar  Cap  Element  35. 
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Figure  3.22x.  Axial  Stress-Spar  Cap  Element  37 
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Figure  3.22y.  Axial  Stress-Spar  Cap  Element  38 


STRESS  VS  PERCENT  LOHD~BE#Jy 


I  bfcJuD 


Expejrin 

lenpal 

•  *  ■ 

— 

— 

Fi 

ni 

te 

El  feme 

nt 

Anial 

/si 

s 

.... 

2 

r 

: 

i 

> 

"2 

r 

; 

; 

s 

> 

s 

s 

> 

s 

'y* 

r 

s 

r 

S' 

”>T'" 
r  \ 

.... 

! 

. . . . 

•  “■ 

PERCENT  LOAD 


Fiqure  3.22z. 


Axial  Stress-Spar  Cap  Element  39 


sgoois'aiSGJOo 

—  ru  .-•-I  T  in  I--  CO 


PERCENT  LOAD 


Figure  3.22aa.  Axial  Stress-Spar  Cap  Element  40 
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Fiqure  3.22dd.  Axial  Stress-Spar  Cap  Element  71 
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Figure  3.22ff.  Axial  Stress~Spar  Cap  Element  73 
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Figure  3. 22gg.  Axial  Stress-Spar  Cap  Element  74 


The  stresses  in  spar  caps  35-40,  and  69-74 
(see  Figure  3.21e)  are  considered  in  Figures  3.22v-gg.  These 
results  do  not  seem  to  have  much  of  a  pattern;  some  of  the 
analytical  vs.  experimental  comparisons  (Bar  elements  35  and 
69,  for  example)  are  quite  close.  Other  comparisons  (Bar 
elements  36  and  70,  for  example)  are  fairly  good  for  lower  load 
levels.  And  still  other  comparisons  are  not  good  at  all. 


RCPUCA  UINC  ✓  MOM-LINEAR  ANAUVSIS  ✓  DAMAGE  IN  CHORD  BAV  1-3  SPAN  BAV 
CASE  3 


UWCR  KIN 


BtPUCA  UINC  /  NON-LINEAR  ANALYSIS  /  DAMAGE  IN  CHORD  BAV  1-3  SPAN  BAV 
CASE  3 


•23b.  Contours  of  Equivalent  Stress  -  40%  Load. 


SECTION  4 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

An  experimental  test  facility  has  been  developed  for 
performing  realistic  simulation  of  flight  loads  on  sections  of 
aircraft  wing  structures.  The  test  facility  is  a  self-contained 
flight  loads  simulation  fixture  consisting  of  a  structure  which 
applies  loads  to  one  end  of  a  test  specimen  and  supports  the 
specimen  at  the  other  end,  a  hydraulic  system  which  imposes 
loads  through  hydraulic  actuators,  and  a  control  system  to 
provide  the  test  operator  with  a  convenient  means  for 
controlling  the  application  of  the  actuator  loads  during  a 
test.  The  test  facility  was  designed  to  operate  within  the 
confines  of  the  Vertical  Gunfire  Facility  at  the  Wright-Patterson 
Air  Force  Base.  The  test  structure  is  self-reacting,  and  thus 
imposes  no  loads  to  the  Vertical  Gunfire  Facility  other  than 
dead  weight.  The  framework  has  been  designed  to  be  open  in  the 
area  of  the  specimen  so  that  the  air  flow  and  projectile  impact 
capabilities  of  the  existing  facility  are  not  interfered  with. 

The  operation  of  the  test  facility  is  performed  through  a 
control  console  located  at  a  site  remote  from  the  physical  test 
area.  The  operator  can  impose  realistic  flight  loads  (bending 
moment,  shear,  torque,  etc.)  to  a  test  specimen  by  adjusting 
dials  which  control  the  load  imposed  on  the  load  frame  by 
hydraulic  actuators.  Proportional  increases  in  the  actuator 
loads  can  be  controlled  conveniently  by  the  operator;  this 
capability  is  useful  for  determining  the  residual  strength  of  a 
damaged  wing  section,  for  example. 

Experiments  were  performed  on  several  replica  test 
specimens,  both  undamaged  and  damaged,  using  the  experimental 
test  facil  ty.  On  the  basis  of  these  tests,  it  is  concluded 
that  the  experimental  test  facility  does  in  fact  provide  the 
survivability/vulnerability  engineer  with  a  convenient  means 
for  imposing  realistic  flight  loads  on  sections  of  aircraft 
structures. 


In  addition  to  the  experimental  test  facility,  an 
analytical  technique  was  developed  for  predicting  the  internal 
load  distributions  of  ballistically  damaged,  multiple  load 
path  aircraft  wing  structures.  The  analytical  phase  consisted 
of  the  modification  of  the  three-dimensional  nonlinear  finite 
element  program  MAGNA  to  include  membrane,  bar,  and  shear  panel 
elements,  the  development  of  a  preprocessor  to  automatically 
generate  the  data  for  the  finite  element  program  given  only  an 
abbreviated  set  of  data,  and  the  development  of  a  postprocessor 
to  display  graphically  the  results  of  the  finite  element 
program.  At  the  suggestion  of  the  Air  Force,  the  analytical 
procedure  was  based  on  the  use  of  simplified  finite  element 
models  of  wing  structures  consisting  of  two-dimensional  membrane 
elements  for  the  skins,  two-dimensional  shear  panel  elements 
for  the  spar  and  rib  webs,  and  one-dimensional  axial  bar 
elements  for  the  spar  and  rib  caps.  Using  these  elementary 
finite  elements,  a  model  of  a  wing  structure  is  built  using  one 
element  per  bay.  The  use  of  the  analytical  procedure  for 
predicting  the  response  of  aircraft  wing  sections  (the  same 
sections  which  are  tested  in  the  experimental  facility)  is 
carried  out  using  three  computer  programs  as  shown  in  Figure 
2.9.  The  programs  are  linked  together  by  the  creation  of 
files;  that  is,  the  output  of  one  program  is  used  as  input  to 
the  succeeding  program. 

The  three  parts  of  the  analysis  procedure  (preprocessor, 
finite  element  program,  and  postprocessor)  all  operate  as  they 
were  intended  to  operate.  The  preprocessor  generates  finite 
element  models  for  undamaged  and  damaged  wing  structures.  This 
is  an  interactive  program  which  requires  only  a  relatively 
small  amount  of  data  defined  in  terms  convenient  for  the  user. 
The  primary  limitation  of  the  preprocessor  is  that  it  is 
restricted  to  generating  finite  element  models  using  only  bar, 
membrane,  and  shear  panel  elements  as  suggested  by  the  Air 
Force.  The  MAGNA  finite  element  program  is  a  comprehensive 
finite  element  program  specially  developed  to  solve  complex 


nonlinear,  static  or  dynamic  structures  problems.  It  has  the 
capability  of  solving  problems  of  a  more  general  nature  than 
those  generated  by  the  preprocessor.  The  postprocessor  is  a 
tool  for  presentation  of  the  output  of  the  finite  element 
program  in  convenient  graphical  form  such  as  undeformed  and 
deformed  geometry  plots,  and  contour  and  relief  plots  of 
displacements,  stresses  and  strains. 

All  of  the  tests  performed  with  the  experimental  test 
facility  were  simulated  with  the  analytical  prediction  tool. 

The  comparison  between  the  experimentally  obtained  and 
analytically  predicted  stresses  in  undamaged  and  slightly 
damaged  replica  test  specimens  were  satisfactory,  especially 
when  the  response  was  in  the  linear  range.  However,  when  a 
substantial  amount  of  damage  was  present  and/or  when  the 
response  was  nonlinear,  the  experimental  and  analytical  results 
differed  substantially.  The  University  of  Dayton  thinks  that 
the  discrepancies  in  the  results  are  due  to  the  use  of 
simplified  bar/membrane/shear  panel  finite  element  models. 

This  type  of  finite  element  model  is  generally  adequate  for 
undamaged  wings  because  they  act  like  box  beams.  However,  the 
local  response  about  a  damaged  area  in  a  wing  model  cannot 
be  adequately  represented  by  such  a  simplified  model. 

On  the  basis  of  the  results  generated  it  can  be  concluded 
that  the  accurate  experimental  and  analytical  prediction  of 
the  stress  distribution  in  damaged  aircraft  structures  is 
feasible.  However,  certain  refinements  to  the  testing  apparatus 
and  to  the  analytical  modelinq  procedure  are  recommended  to 
improve  the  correlation  between  experimentally  and  analytically 
obtained  data.  The  primary  area  for  improvement  to  the  testing 
apparatus  is  the  method  for  attaching  test  SDecimens  to  the 
loading  frame  and  to  the  reaction  fixture.  It  is  recommended 
that  the  attachment  brackets  be  redesigned  to  prevent  relative 
motion  between  the  specimen,  and  the  reaction  and  loading  frames. 
This  design  change  is  necessary  for  improving  the  accuracy  of  the 
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displacements  which  are  obtained  experimentally.  Every  attempt 
should  be  made  to  measure  accurate  displacements,  since 
displacements  are  the  most  accurate  data  predicted  analytically 
by  the  finite  element  approach.  It  is  also  recommended  that 
additional  strains  be  recorded  to  obtain  a  more  complete 
description  of  the  stress  state  in  a  damaged  test  item,  for 
comparison  to  analytically  predicted  data. 

The  primary  area  for  improvements  of  the  analytical 
prediction  of  the  response  of  damaged  wings  is  the  method  for 
finite  element  modeling.  While  the  use  of  bars,  membranes, 
and  shear  panels  are  certainly  sufficient  for  representing  the 
behavior  of  undamaged  wings  which  respond  essentially  as  box 
beams,  the  results  obtained  indicate  that  this  approach  does 
not  produce  accurately  the  stress  response  of  damaged  wings. 
Thus,  it  is  recommended  that  a  study  be  conducted  to  determine 
alternative  modeling  procedures  for  predicting  damaged  wing 
response. 


REFERENCES 


Roth,  F.  J.  ,  R.  J.  Dominic,  and  D.  Bauer,  "Structural 
Flight  Loads  Simulation  Capability  -  Interim  Technical 
Report  -  Load  Fixture  Design,"  UDRI-TR-77-40 ,  University 
of  Dayton,  June  1977. 

Bruner,  T.  S. ,  M.  P.  Bouchard,  M.  J.  Hecht,  and  F.  K. 
Bogner,  "Structural  Flight  Loads  Simulation  Capability  - 
Structural  Analysis  Computer  Program  User's  Manual," 
UDR-TR-80-73 ,  University  of  Dayton,  June  1980. 

Fiscus,  I.  B. ,  "Replica  Test  Specimen  Design,"  UDR-TR- 
80-25,  University  of  Dayton,  January  1980. 

Brockman,  R.  A.,  "MAGNA  Computer  Program  User's  Manual," 
UDR-TR-80-107 ,  University  of  Dayton,  November,  1980. 


176 


■*U.S. Government  Printing  Office:  1981  —  757-002/377 


